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It is by now well known that various condensed matter systems may be used to mimic many of the 
kinematic aspects of general relativity, and in particular of curved-spacetime quantum field theory. 
In this essay we will take a look at what would be needed to mimic a cosmological spacetime — to 
be precise a spatially flat FRW cosmology — in one of these analogue models. In order to do this 
one needs to build and control suitable time dependent systems. We discuss here two quite different 
ways to achieve this goal. One might rely on an explosion, physically mimicking the big bang by 
an outflow of whatever medium is being used to carry the excitations of the analogue model, but 
this idea appears to encounter dynamical problems in practice. More subtly, one can avoid the 
need for any actual physical motion (and avoid the dynamical problems) by instead adjusting the 
, propagation speed of the excitations of the analogue model. We shall focus on this more promising 

route and discuss its practicality. 
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I. MOTIVATION 



O 

in 

It is by now well known that various condensed matter systems may be used to mimic many of the kinematic aspects 
of general relativity, and in particular of curved-spacetime quantum field theory 0, 0, U 13 • Prior work has focussed 
largely on black holes and event horizons, but there are also interesting cosmological issues that can be addressed. In 
this essay we will take a look at what would be needed to mimic a cosmological spacetime — to be precise a spatially 

CJ^ flat FRW cosmology — in one of these analogue models. There are two quite different ways of trying to achieve the 

5^ \ same goal: 

• One might rely on an explosion, physically mimicking the bigbang by an outflow of whatever medium is being 
. used to carry the excitations of the analogue model (see e.g. |3)- Unfortunately, this idea appears to encounter 

dynamical problems in practice, problems not inherent to the type of geometry being reproduced. 

• More subtly, one can avoid the need for any actual physical motion (and avoid the dynamical problems) by 



instead adjusting the propagation speed of the excitations in the analogue model. 
The physical metric we will be trying to emulate is that of a spatially flat FRW cosmology 

d4 RW = -c 2 di 2 + a{tf dx 2 , (1) 

where a(t) is the scale factor of the universe as a function of time. In contrast the analogue models generically provide 
effective metrics of the form El Q 

dative = ^jA { ~ [cl (*, x) - v 2 (t,x)] dt 2 ~2v(t,x)- dt dx + dx 2 } . (2) 

Here v(t,x) is the physical velocity of the medium, c s (t,x) is the propagation speed of whatever excitations we are 
interested in studying, and the conformal prefactor p(t,x)/ c s (t,x) depends both on the dimensionality of spacetime 
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and [to some extent] on the specific choice of the analogue model. In the form presented above the conformal factor 
is appropriate to either ordinary sound in a classical fluid Q, 0, 0, 0], or to phonons in a BEC 6], both in (3+1) 
dimensions. 

We now want to consider the fundamental question — given the availability of "effective metrics" of type J2J), how 
close can one get to reproducing a FRW metric of type (JIJ . That this question is non- vacuous can be deduced from the 
observation that in the black hole context "effective metrics" of type J2J never exactly reproduce the Schwarzschild 
geometry ||, they can do so only up to an overall conformal factor. We wish to check, in particular, if the same 
phenomenon shows up in a cosmological context. 

II. EXPLOSIONS 

As mentioned in the introduction, one plausible approach is to rely on an actual physical explosion in the medium to 
mimic the big bang [jj. Start with the FRW metric Q and substitute 

z = b(t)x; dz = b(t) dx + bit) x dt; b(t) dx = dz- zdt; (3) 

b(t) 



Then, introducing a Hubble- like parameter, 



*(«) = (4, 



the FRW metric ^ is transformed to 



as FRW 



c 2 -^ H b (tf ? 



,2 



dt 2 - 2 % H b (t) z -dz dt + % dz 2 . (5) 
b z b 



This metric is equivalent to the required "effective metric" form [equation PJl]. provided we take 

_, . . _ b p a 2 , . 

v(t, z) tw) H b (t) z; c s t*t - C ; — w» (6) 

a c s b z 

The continuity equation for the medium, 

p + V-( P v) = 0, (7) 
particularized to the flow field v(t, z) H b {t) z implies 

p + 3pH b (t)=0; p (t)ocp^. (8) 

Using the last equivalence in Eq. © we can now fix the behaviour of c s (t) as a function of a(t) and b(t) 

Equations ijHJ and 0, together with the first equivalence in Eq. ©, completely fix the relation between the hydro- 
dynamical quantities (p(t), v(t, z), c s (t)) and the cosmological solution parameters (a(t), b(t), c). 

Although the "explosion route" just discussed seems promising, one should note that because of the linearly rising 
velocity field v = H b z, this particular realization of a FRW effective geometry is guaranteed to possess an apparent 
horizon, a spherical surface in which the speed of the fluid surpasses the speed of sound. From a dynamical point of 
view, this might introduce many practical problems not intrinsically inherent to the type of geometries we are trying 
to reproduce. In this sense, we view the use of an exploding medium as not being a particularly useful analogue for 
an expanding FRW universe. 
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III. VARYING PROPAGATION SPEED 

A much better analogue of an expanding FRW universe can be obtained by keeping the medium at rest and instead 
varying the propagation speed of the excitations. In this case the continuity equation implies that p{t, x) is a constant, 
so that we can rescale the metric by a constant factor cq/ p (where Cq is any convenient reference speed) obtaining 

d S 2 ffcctivc = -^L {- c l(t) dt 2 + dx 2 } = - Co c s (t) dt 2 + -BL dx 2 . (10) 



Now introduce a pseudo-time r, related to laboratory time via dr = dt y/c s (i)/co. Then 



ds rffoctivc = - c o dr2 + ^) df2 - ( U ) 
Consequently we completely reproduce a FRW cosmology provided we identify 

c dr » c dt FRW ; — a(r) 2 . (12) 

That is, an expanding universe corresponds in the effective geometry to a decreasing propagation speed. A tricky point 
is that c s is still presented in terms of laboratory time t as the speed ||da?/d£||. However in terms of the pseudo-time 
r the excitations propagate at speed 



c s (t) = 



dx 


d* 


dx 


d7 


~ d7 


dJ 



c s (i) 



c,(t) = Vco c a (t). (13) 



Having obtained the general analogue of the flat FRW metric it is interesting to investigate the analogue equivalent 
of the inflationary solution. We can start by re-writing the analogue Hubble factor as 

a(Tj 2 y Co dr 2 y c s c s 

where the prime represents derivative with respect to the pseudo-time and the dot with respect to the laboratory 
time. The inflationary solution is then easily obtained by using a exponential law in pseudo-time r 

a(T)=e HT ; a(t) = H t; (15) 

corresponding to a power law in physical time for the speed of sound (measured in physical time) 

=■(') -7%i-wr < 16 ' 

In summary the nice feature of the variable propagation speed route to a FRW analogue is that it is relatively clean; 
there is no moving medium which might impact the physical boundaries of any experimental apparatus and one can 
instead focus on what we feel is the central issue — how to physically manipulate the propagation speed c s (or its 
pseudo-time equivalent c s ). 

IV. SUITABLE PHYSICAL MECHANISM: FESHBACH RESONANCE IN A BEC 

From the preceding discussion it is clear that one attractive route to simulating a FRW cosmology in an effective 
geometry is by rapid manipulation of the propagation speed of whatever excitations we might wish to focus attention 
on. But how is that to be accomplished? There is one particular medium, currently the center of much experimental 
interest, for which an appropriate mechanism has been demonstrated to exist. Here we are referring tho the use of 
Feshbach resonances in Bose-Einstein condensates [BECs] . 

BECs are promising analogue systems of gravity. They have been extensively studied 0, 0, S @ m relation to the 
possibility of simulating black hole geometries and Hawking radiation. The basic equation used in describing the 
condensate is the Gross-Pitaeveskii [GP] equation 

d { Ti? \ 4 /j 2 

-m— V^,x) = (-— V 2 + 14xt(x)+ K |^(i,x)| 2 W(t,x); k=^- (17) 
at \ 2m I m 
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Here i/)(x, t) is the (classical) wave function of the condensate, (\ip(x,t)\ 2 = n =thc particle density of the condensate), 
V^xt is the trapping potential and a is the s-wave scattering length for the atoms (which have mass m). [Actually it 
is conventional to use the symbol a for the scattering length of a BEC condensate. We adopt this unusual notation to 
avoid any confusion with the scale factor of the FRW metric] The generic analogue model can be easily be obtained 
by considering the propagation of excitations in the condensate. These arc described by a wave function which closely 
resembles that of a scalar field in a curved spacetime described by the metric (0 with p replaced by n m and where 
c 2 s = nn/ra (see for a detailed derivation). The key point is that the propagation speed is proportional to the 
scattering length 

c 2 s oc a. (18) 

Let us now qualitatively explain we can simulate cosmological expansion within this model. We start with a condensate 
in a stationary state described by a constant background density in a sufficiently large volume. This is a solution of 
the Gross-Pitaevskii equation. For this, one needs to have a potential that reproduces a sufficiently large hard-walled 
box. Using a Feshbach resonance it is now possible to change at will the scattering length a in the condensate. If one 
now decreases the value of the scattering length in a sufficiently slow manner (the timescale will be discussed below) , 
then at the same time the value of the speed of sound will decrease. The analysis of the previous Section then shows 
that this leads to an effective FRW geometry Fluctuations of the condensate will propagate in an effective metric 
which is an analogue to a spatially flat FRW geometry. 

A subtle point is that when the time dependence is introduced in the system by a time- varying scattering length one 
must be sure to work in a regime where the background configuration is "instantaneously" reacting to theses changes. 
Thus we must assure that the GP equation (from which the analogue gravity framework is derived) holds at each 
instant of time. The validity of the GP in describing the Bose-Einstein condensate is related to the validly of several 
crucial assumptions, generally stated to be the "mean-field" approximation and the dilute gas approximation. It is 
nevertheless important to note that in a dynamical situation a third approximation, which we can call "Markovian" 
approximation, is implicitly assumed. 

The Markovian approximation is related to the fact that in dynamical situations the two-body time-dependent scat- 
tering matrix can have a complicated form due to the "memory" of the system. In these situations the system is 
described by a GP-like equation where the interaction term includes a "delay term" described by an integration over 
time. The necessary assumption in order to have a Markovian description of the dynamics (which together with the 
two previous approximations leads to the GP equation) is then that the timescales on which external parameters are 
changing are longer that the two-body collisional duration. That is, longer than the timescale over which a single in- 
teraction happens. (Reduced to the bare bones, we are asking that the scattering length does not change significantly 
during the period when a pair of atoms are interacting.) 

We can estimate the two-body collision time by a simple calculation. All we need is the typical size of the region of 
strong interaction of two atoms in the condensate, and the typical speed with which they move. The first quantity 
can be assumed to be of the order of the Van der Waals length. This length is basically the size of the region of strong 
interaction: For r <C X v dw the scattering wave function oscillates rapidly due to the strong interaction potential. In 
alkali ground state interactions, X v dw is typically of order of few nanometers, X vc iw ~ 1 nm - 

Regarding the typical speed of the atoms, this is set by the de Broglie momentum generated by the trap confinement: 
p = h/R and v — p/m. We shall assume a trap of typical size of ten microns. For a typical Bose-Einstein condensate 
one then gets an interaction timescale 

u = W = X vdw mR w 10 _ 6 (19) 
v h 

so it would appear that a microsecond is the shortest timescale allowed for the change in a. 



V. ANALOGUE COSMOLOGICAL PARTICLE CREATION 



Now that the above analysis has shown that there is a regime for which a varying scattering length can be used to 
simulate a FRW-like effective metric one may wonder about the behaviour of the fluctuations on such a time-dependent 
background. 

The equation satisfied by the quantum fluctuations is, in the acoustic approximation (that is for long wavelengths), 
that of a massless scalar field over an expanding background and therefore, it will lead to cosmological pair production 
of particles. Interestingly it was recently demonstrated 01 that not all analogue models of gravity are suitable for 
simulating particle creation from the quantum vacuum. Indeed, it may happen that even if the classical equation 
for the perturbations on the background resembles that for a field in curved spacetime, nonetheless one might fail to 
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mimic quantum particle creation because the commutation relations of the analogue field are not the correct one for 
the physical quantum field. In this regard it is useful to note that it can be explicitly shown [lOj that in the case of 
BECs the structure of the commutator is the correct one. 

We then expect to be able to produce quantum excitations with several frequencies, typically determined by the 
rapidity of the change in the scattering length. In particular if i m ; n is the shortest timescale over which we can 
physically drive the system then ^ p0 ak ~ 1/imin is the peak frequency of the created quasi-particles spectrum. It 
is also important to check that this frequency corresponds to a wavelength shorter than the physical size of the 
condensate R. An easy way to do this is to confront the interaction time ti given in Eq. I|19|) with the crossing time 

R/c. 



tsize h £ 

where we have introduced the healing length of the condensate, defined to be £ = h/(mc s ). For typical BEC system 
£ 0.1-10/xm (assuming an average value of the scattering length of a few nanometers) so 

A. » l(T 2 -l(n 4 . (21) 

^-sizc 

Thus there is a viable window of timescales (for the time dependence of the scattering length) for which both the GP 
equation holds, and the quasi-particles produced have wavelengths shorter than the physical size of the condensate. 
The possibility of simulating cosmological particle creation is a significant advance that deserves further investigation, 
and might be remarkably important in the future. The simulation of inflationary scenarios, and the relevant particle 
creation, could lead to a better understanding of the generation of primordial inflationary perturbations and their 
role in the generation of large scale structure. Moreover it should be noted that the quasi-particles generated in this 
way will be characterized by a phononic dispersion relation lu 2 = c 2 k 2 only at long wavelengths with respect to the 
healing length of the condensate. In general for wavelengths comparable with the healing length of the condensate 
the Bogoliubov dispersion relation 



_^ fc2 + (A fc2 ) , (22) 

will hold • This is particularly interesting due to the recent intense debate about the detectability of inflationary 
spectrum due to these kind of dispersion relations 0, . 



VI. CONCLUSIONS AND PROSPECTS 



Our major conclusions are four-fold: 

• At a theoretical level, it is clear that any mechanism for changing the propagation speed in a stationary medium 
is from a mathematical perspective equivalent to working in an expanding (or collapsing) spatially flat FRW 
universe. 

• At an experimental level, the use of a Feshbach resonance in a Bose-Einstein condensate yields a way of influencing 
the scattering length, and hence the propagation speed, without changes in condensate density — and more 
importantly without violating the approximations made in deriving the Gross-Pitaevskii equation on which the 
entire "effective metric" approach to phonon propagation in BECs is based. 

• At a very practical engineering level, the relevant parameters seem to be well within our technological horizon. 

• Building an analogue FRW cosmology (suitable for testing semiclassical quantum effects) seems considerably less 
problematic than building an analogue black hole 0,0; Ell- 
in summary: The prospects for direct laboratory simulation of an expanding universe, and consequent cosmological 
particle production are very good. Given the relatively small number of experimental tests of curved space quantum 
field theory, any progress along these lines is important. Our long range goal is to turn at least some aspects of 
cosmology into a laboratory science, not just an observational science. 
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